The nonlinear elasticity of an a-helical polypeptide 
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We study a minimal extension of the worm-like chain to describe polypeptides having alpha- 
helical secondary structure. In this model presence/absence of secondary structure enters as a 
scalar variable that controls the local chain bending modulus. Using this model we compute the 
extensional compliance of an alpha-helix under tensile stress, the bending compliance of the molecule 
under externally imposed torques, and the nonlinear interaction of such torques and forces on the 
molecule. We find that, due to coupling of the "internal" secondary structure variables to the 
conformational degrees of freedom of the polymer, the molecule has a highly nonlinear response to 
applied stress and force couples. In particular we demonstrate a sharp lengthening transition under 
applied force and a buckling transition under applied torque. Finally, we speculate that the inherent 
bistability of the molecule may underlie protein conformational change in vivo. 

PACS numbers: 87.10.+e,87.14.Ee,82.35.Lr 
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I. INTRODUCTION 

The study of the mechanical properties of individ- 
ual biopolymers serves as an important laboratory to 
probe polymer physics at the length scale of a single 
chain and further elucidates the biological processes in 
which these molecules take part 0,0- For example the 
biologically fundamental processes of DNA replication, 
transcription, the regulation of transcription rely on the 
DNA protein interactions involving the mechanical defor- 
mation and microstructural modification of DNA both 
at long length scales as well as at the scale of indi- 
vidual base pairs. Recent advances in the experimen- 
tal manipulation of individual biological macromolecules 
has opened a new window on these processes and al- 
lows for the direct quantitative test our understanding 
of the mechanical properties of these macromolecules in 
thermal equilibrium. These single-molecule manipula- 
tion experiments haveprobed the mechanical properties 
of not only DNA@, S II H H 0, H but also a vari- 
ety of biologically important macromolecules including 
polysaccharridesKl llCj. and giant proteins such as titin 
|llllT^ . IT^| and tenascin [TJ. The better theoretical un- 
derstanding of protein mechanics will enhance the inter- 
pretation of protein force spectroscopy, which may even 
shed light on protein folding pathways [l5j although this 
latter point appears to be somewhat controversial [16l 
Regardless, understanding the mechanical properties of 
proteins is fundamental to elucidating the allosteric or 
conformation changes that many proteins undergo as 
part of their biological activitv flSl fl9t 12(1 l2l| . 

These last examples demonstrate the feasibility of the 
direct mechanical manipulation of single proteins. How- 
ever, modelling the mechanical properties of these atyp- 
ically large proteins or, for that matter, any entire pro- 
tein is a daunting task since such molecules have complex 
structures that result from a number of local and non- 
local interactions along the polymer chain. In order to 
make quantitative progress in the interpretation of these 



protein manipulation experiments it appears to be useful 
to first understand in more detail the mechanical proper- 
ties of simpler polypeptide-based structures. A natural 
candidate for such a simpler structure is a protein sub- 
domain of one secondary structure. Here too there is 
experimental input: single molecule force spectroscopy 
via AFM has been used to directly probe alpha-helical 
polypeptides |22| as well as synthetic polymer chains with 
a local helical structure such as PEG 
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In order to understand the mechanical properties of 
proteins in general and alpha-helical polypeptides in par- 
ticular it is necessary to develop a new minimalistic 
model that incorporates both the conformational fluc- 
tuations of the polymer backbone and localized struc- 
tural transitions of the constituent monomers. In other 
words it is necessary to augment simple models of the 
statistical mechanics of the peptide backbone with new 
terms to account for the presence of secondary structure 
along the chain. Furthermore, we must allow for interac- 
tion between the degree of local secondary structure and 
the conformational degrees of freedom of the polypetide 
backbone. 

In this paper we examine the predicted mechanical 
properties of such a minimalistic model of an alpha- 
helical polypeptide in which we allow the interaction of 
the local secondary structure of the chain with its con- 
formational degrees of freedom. We treat the local pres- 
ence or absence of secondary structure as a two-state 
(Ising-like) variable along the chain backbone, which is 
itself described by a set of local tangent vectors to the 
chain. In order to make this simplification, we coarse- 
grain the polymer so that each independent monomer can 
be unambiguously assigned a state of secondary struc- 
ture. This requires us to consider a model comprised of 
coarse-grained chain segments (i.e. monomers) each con- 
sisting of about three amino acids. The interaction be- 
tween the internal, secondary structure variables and the 
conformation of the polymer chain as described by the 
set of backbone tangent vectors is effected by the pres- 
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ence of a bending modulus of the backbone whose value 
depends on the local state of secondary structure. When 
these three amino acids making up one model monomer 
adopt a local configuration consistent with alpha-helical 
secondary structure, the hydrogen bonding between these 
amino acids [24| renders that segment of the chain signif- 
icantly stiffer than the same polymer without the locally 
ordered secondary structure. Thus the bending energy 
associated with the local change in the backbone tangent 
vectors is higher in regions having alpha-helical structure 
than in regions locally adopting a random coil configu- 
ration. Similar models can and have been applied to 
study the mechanical properties of DNA and have been 
discussed in the current context as well[25l I2H l27j. 

To qualitatively characterize our results presented be- 
low, we note that, due to the presence of the internal 
state variables representing secondary structure along the 
chain and their control over the local chain bending mod- 
ulus, the alpha- helix is predicted to have a highly non- 
linear response to both bending torques and to exten- 
sional forces. Under small externally applied torques, 
the molecule will deform so that its thermally averaged 
chain contour takes the form of the arc of a circle and 
the torque necessary for bending the molecule through a 
given angle grows linearly with that angle. The molecule 
deforms roughly as a flexible rod. At a critical torque, 
however, the secondary structure of the molecule is lo- 
cally disrupted producing a small length of the backbone 
with a much softer bending modulus. The total curva- 
ture that had been uniformly distributed along the back- 
bone becomes localized in the anomalously soft region 
produced by the disruption of the secondary structure 
and the torque required to enforce the curvature of the 
molecule drops precipitously. 

The long range goals of this sort of modelling go be- 
yond the interpretation of the emerging experiments on 
the mechanical properties of polypeptides having alpha- 
helical secondary structure. By understanding the me- 
chanical properties of the constituent elements of a pro- 
tein it should be possible to develop a lower dimen- 
sional representation of protein mechanics. In place of 
the atomic coordinates of the backbone carbons and the 
positions of the various amino acid residues, one may de- 
scribe protein domains (having definite secondary struc- 
ture) as a space curve having some nonlinear extensional 
and bending compliances that may be computed in terms 
of a few energy scales determined either from experiment 
or simulation. Using three-dimensional protein struc- 
tural data and such a nonlinear elastic model for each 
structural element of definite secondary structure, one 
can attempt to build mechanical models of entire pro- 
teins that, due to their highly reduced number of degrees 
of freedom are more tractable for numerical investigation 
than those based on all atom simulations. 

From the study of those models one may be able to 
extract low energy conformational pathways and thus 
make predictions regarding protein allostery. For ex- 
ample, from a combination of native-state protein struc- 



tural data and the calculated nonlinear elastic properties 
of alpha-helical protein domains it may thus be possi- 
ble to predict the mechanical properties of the alpha- 
helical coiled coil region in myosin II |2S| or alpha-helix 
rich proteins such as spectrin as probed by mechanical 
unfolding experiments [2^. l3fj| . Further data on pro- 
tein conformational change is available from numerical 
simulation |3ll |32| . 

The remainder of the paper is organized as follows. 
In section [H] we introduce the alpha-helix Hamiltonian 
based on a combination of the worm-like chain and the 
helix/coil model. Using this model we calculate the re- 
sponse of the chain to bending torques in section UlI Al 
We then take up the problem of the extensional compli- 
ance of the alpha-helix in section Till Bl before summariz- 
ing the results and discussing possible experimental tests 
of the theory in section llVl 

II. THE HELIX-COIL WORM-LIKE CHAIN 

The worm-like chain (WLC)|U|3| is the fundamen- 
tal coarse-grained model for a polymer at length scales 
shorter than its thermal persistence length. This model 
describes the single-chain polymer statistics in terms of 
a quadratic Hamiltonian that associates an energy cost 
with chain curvature by introducing a bending modulus 
k. In terms of a discretized chain model described by the 
set of monomeric tangent vectors tj, i = 0, . . . , N— 1 with 
N the degree of polymerization, the WLC Hamiltonian 
may be written as 

N-l 

H W LC = K [1 - (*< ' *Vl)] ■ (1) 

i=0 

The effect of this bending energy is to enhance the statis- 
tical weight of straight chain configurations on a length 
scale of k monomers equal to the thermal persistence 
length. Here and throughout this paper we take k^T = 1. 
It may be easily checked that this length is equal to the 
arc length (measured in monomer lengths 7) of the poly- 
mer chain over which the chain tangent vectors thermally 
decorrelate. At length scales much longer than kj the ef- 
fect of this bending energy is minimal and the equilibrium 
statistics of the polymer become controlled by a combi- 
nation of intrachain collisions and chain configurational 
entropy |35j with a renormalized Kuhn length. In theta 
solvent one finds in the limit of very long worm-like chains 
the radius of gyration to be given by (iS) = 2kj 2 N 
where 7 is the monomeric length. 

The single-chain response to externally applied ten- 
sion F has also been exhaustively researched within the 
WLC description of the polymer [36|, H3, |3g • The fun- 
damental result of this work is that one may determine 
the extensional compliance of the molecule as a function 
of applied force. This compliance is defined as d(L)/dF, 
the derivative of the equilibrium chain length (L) with re- 
spect to the applied force F. At low forces it is essentially 
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FIG. 1: Schematic figure of an alpha-helical polypeptide and 
its schematic representation in terms of the Ising-like sec- 
ondary structure variables (open circles for random coil seg- 
ments and filled ones for alpha-helical ones) and the tangent 
vectors to the segments of the chain (denoted by arrows). 



constant reflecting the standard, quadratic reduction in 
chain configurational entropy associated with long, flex- 
ible polymers. At high forces, however, the compliance 
goes to zero as i* 1-3 / 2 due to the fact that the chain has 
a finite length at even arbitrarily high forces. The char- 
acteristic form of the approach of the compliance to zero 
in the high force limit is controlled by the pulling out of 
small, transverse thermal fluctuations of the chain and 
thereby recovering the arc length stored in them to in- 
crease (L). 

While the WLC is a highly successful model to describe 
the force extension properties of a number of biopolymers 
such as DNA, it is clear that it is not sufficient to properly 
describe these molecules under large enough tensions. At 
larger tensile stresses, details of the internal structure of 
the molecule become important for the understanding of 
conformational properties of the molecule. For example, 
under large enough stresses the double-helix structure of 
DNA (B-DNA) can be unwound allowing each monomer 
to lengthen by a factor of about 1.85 |39|. To account 
for such (two-state) internal degrees of freedom along 
the chain, workers have employed the helix/coil (HC) 
model 4X| . This model has been used to study a class of 
protein conformational transitions [4 ll |42| in solution and 
under tension|43j. 

The HC model Hamiltonian, which is used to model 
these structural transitions can be reduced to it simplest 
form by assuming that the local structure of the chain can 
be described by a set of two-state variables Sj = ±1, % = 
0, . . . , N. For the alpha-helical chains of current interest 
we regard these two states as the local conformation of 
the monomer in its native, alpha-helical state (s = +1) 
and in a disordered, random coil state (s = —1). The 
statistics of this set of two-state variable is controlled by 
the Hamiltonian: 



N-l 



N 

H HC = e w /2 - s iSi+1 ) - h/2j2(** - !)■ (2) 

i=0 i=0 

It is immediately clear that above Hamiltonian also de- 
scribes a one dimensional, ferromagnetic Ising chain. The 



interpretation here is somewhat different. The energy h 
playing the role of an external magnetic field in the Ising 
system now represents the free energy cost per monomer 
to be in the non-native (i.e. random-coil) state. This 
term is thus controlled by a combination of the chem- 
istry of the monomeric residues and solvent quality; its 
calculation from fundamental solution chemistry is be- 
yond the scope of the current work, however we will at- 
tempt to estimate its magnitude based on experiment. 
Clearly this constant is at least of order unity since the 
protein domain under investigation is assumed to have an 
alpha-helical secondary structure in thermal equilibrium. 

The first term in Eq. [5] plays the role of the nearest 
neighbor ferromagnetic coupling in the Ising interpreta- 
tion of the Hamiltonian. In its current interpretation, 
e w is free energy cost of a domain wall in the sequence 
of helix (s = +1 ) and random coil (s = —1) sites. In 
the helix/coil literature it is also referred to as the nat- 
ural logarithm of the "chain-cooperativity" parameter. 
By adopting a native state configuration, a monomer 
presents hydrogen bonding sites to its neighbors. If those 
neighbors are also in their native state, these hydrogen 
bonds further lower the free energy of the system via this 
nearest neighbor cooperative effect. If, however, one of 
the neighboring monomers of an alpha helical monomer 
is in its random-coil state, such hydrogen bonding is not 
possible and the total free energy of this domain wall 
configuration is larger than simply the free energy cost 
for one monomer being in the non-native configuration, 
h. 

Finally, we note that since there are multiple hydrogen 
bonds per turn of the alpha helix, it might be reasonable 
to describe the local secondary structure by a g-state dis- 
crete variable where q > 2. Such an analysis changes 
Eq. |21 into an equally tractable one-dimensional g-state 
potts model, but introduces additional unknown param- 
eters. It is thus inconsistent with our goal of exploring a 
minimal model that incorporates secondary structure. 

The coupling of the secondary structure variables to 
the WLC tangent vectors is affected by introducing a 
bending stiffness in the WLC Hamiltonian that depends 
on the local degree of secondary structure. We choose 



k(s) 



k> if s = +1 
k < if s = — 1 



(3) 



Due to the hydrogen bonding between turns of the alpha- 
helix, it is reasonable to expect that k > the bending 
modulus in the native state is significantly larger than 
k < the bending modulus of the chain in the non-native, 
disordered state. We return to the question of determin- 
ing physically reasonable estimates of these quantities in 
the conclusions. By introducing Eq.[3]and combining the 
Hamiltonians in Eas. lll2l we write the full Hamiltonian for 
the coupled system of secondary structure variables and 
chain tangent vectors as the helix/coil worm-like chain 
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Hamiltonian (HCWLC): 



A. Bending 



N-l 



N 



H = e w /2 (! - s «^+i) - V2^(s* - 1) + 

i=0 i=0 

E [1 - (*i ' • (4) 



We note that the system described by Eq. 0] may be 
looked at as two intercalated Heisenberg and Ising 
(s,) magnetic systems. The nearest neighbor coupling of 
the Heisenberg (chain tangent) vector , however, depends 
on the value of the Ising (secondary structure) variable 
between them. A pictorial representation of the system 
along these lines is shown in figure 2] The full Hamilto- 
nian given by Eq. 01 has four constants with dimensions 
of energy: k>, k<, h, e w that can be fit from experiment. 
We return to this point in our conclusions. Finally, we 
note that we have disregarded the twist degree of free- 
dom of the molecule. Such twist degrees of freedom and 
the coupling of twisting and stretching modes of these 
chiral molecules have been explored particularl y w ith re- 
gard to the mechanical properties of DNA [44L l45j . This 
extension of the basic model will be explored in future 
work. 



III. MECHANICAL PROPERTIES 



We now turn to the exploration of the mechanical prop- 
erties of the polymer described by the HCWLC Hamilto- 
nian. We study this problem in different ways. First, we 
consider the response of the molecule to externally ap- 
plied torques by examining the torque required in ther- 
mal equilibrium to enforce a given angular deviation be- 
tween the first and last chain tangent vectors. Second, 
we study the force-extension relations for this molecule 
by calculating the projection of the chains's mean end- 
to-end distance along the direction of an applied force as 
a function of the magnitude of that force. We study this 
extensional compliance for two different cases. In one 
case we assume that the end tangents of the chain re- 
main unconstrained. In the second, we explore the effect 
of applied torque on the extensional compliance of the 
molecule by first constraining the end tangents vectors 
and then applying an extensional force. All these calcu- 
lations are conceivable as individual, single molecule ex- 
periments through the use of e.g. optical and magnetic 
traps. By understanding the dependence of the exten- 
sional compliance on the curvature of the molecule, one 
may gain insight into the mechanical properties of alpha- 
helical domains of proteins that are similarly constrained 
in the protein's native state. 



To consider the bending response of the chain to ap- 
plied torques in thermal equilibrium we first express the 
restricted partition function of the system subject to the 
constraint that the chain tangent deflects by a fixed angle 
ip over its total arc length. It is reasonable to suppose 
that the chain bends in the plane defined by first and 
last chain tangent that are being constrained; to sim- 
plify the calculation we assume that we may examine the 
problem in this two-dimensional subspace. In that case 
the chain tangents are each equivalent to single angle: 
£j_l_l ■ fj — > cos(#i + i — 9i) and we write the restricted 
partition function as 



N 



N 



ZW = ]J E r[{de i e- H 8(e )5{e N -i>), (5) 



i=0 



where the delta functions enforce the constraints on the 
initial and final chain tangents. The Hamiltonian appear- 
ing above is given by Eq. 0] and is a function of all the 
secondary structure variables and backbone tangents. 

To evaluate the partition function it is useful to ob- 
serve the formal equivalence of Eq. [S] to the imaginary 
time propagator of a quantum particle on the unit circle 
|46j |. The partition function above is an imaginary-time 
sliced path integral representation of the transition am- 
plitude for the particle to start at angle 9q = and end at 
angle 9n = if) in N time slices. The imaginary time evo- 
lution operator is simply the exponentiated Hamiltonian 
appearing in Eq. [S] The quantum analogy is somewhat 
complicated by the presence of the secondary structure 
variables; for the current problem the fictitious quantum 
particle has a two-level "internal" variable similar to the 
spin states of a spin 1/2 particle. The state of this par- 
ticle in the angle representation take the form: |s, 6) and 
Eq. [S] can be recast in the form 



E 

S0,SjV = ±l 



Vs N (so,do = 0\T N \s N ,9 N =il>), (6) 



where T is the single-step imaginary time evolution op- 
erator and rj s = e~ h & s ,-i is a factor needed to correct 
the statistical weight of finding the last chain monomer 
in the disordered, non-native state. The remaining sum 
is over the starting and ending secondary structure of the 
chain. 

To make progress it is useful to work in terms of the 
integral angular momentum variables rrii conjugate to the 
angles 9i. In this angular momentum representation of 
the problem we may expand the angle eigenstates as 



\s,6) = \s) 



1 



E 



m— — oo 



-ijn6 



(7) 



In this momentum representation the time evolution op- 
erator is diagonal, i.e. it connects states of the same m 
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only: 



- K> I m [K>] 
- h - K <I m [K<] 



(ms\T\m s ) 

e- h - K <I m [ K< ] 



(8) 



To obtain the above result we have used an identity relat- 
ing the exponentiated cosine to a sum of modified Bessel 
functions of integer order ^3 : 



o Jcos(0) _ 



m— — oo 



Im(J)t 



(9) 



The action of the transfer matrix or imaginary-time evo- 
lution operator can be further simplified by diagonalizing 
it in the remaining 2x2 subspace of secondary structure. 

One may note that the above matrix (Eq. |HJ is non- 
Hermitian reflecting the lack of time-reversal symmetry 
of the underlying Hamiltonian. This absence of time- 
reversal symmetry occurs because the local bending mod- 
ulus between the i th and (i + l) th chain tangents depends 
only on s i; the secondary structure variable to the right 
of the first tangent vector. The absence of microscopic 
left/right symmetry along the chain results in the non- 
Hermitian character of Eq. [S] Minor modifications of 
Eq. 01 generate Hermitian Hamiltonians expressing the 
same physics at length scales larger than the monomer 
size, but we do not pursue such related problems here. 

The transfer matrix Eq. |S] can be diagonalized in the 
space of secondary structure by a similarity transform us- 
ing the matrix U{m) (defined in appendix^); the eigen- 
values of the transfer matrix are 



Ai, 2 (m) 



UJ, 



[f + z m ±y/{l~z m ) 2 + {3 2 z m ], (10) 



where oj m (n) = exp(— K)I m {K) is the transfer matrix el- 
ement for an ordinary WLC and z m = e^ ""^] is 
the ratio of the fugacities of a random coil segment and 
an a-helical segment for a given angular momentum m. 
The quantity (3 appearing in the above expression is 
the exponentiated free energy cost of introducing a do- 
main wall in the secondary structure, helix-coil variables. 
(3 — exp[log[(2) — e w ] with the first and second terms aris- 
ing from respectively the entropic gain and enthalpic cost 
of the creation of a domain wall. 

It is important to note that the partition function of 
the chain may not be reduced simply to the product of 
eigenvalues since doing so presupposes periodic bound- 
ary conditions. While in the thermodynamic limit of 
long chains [N — ► oo) it is indeed permissible to choose 
those nonphysical boundary conditions, if we wish to 
study the dependence of mechanical properties on the 
degree of polymerization it is essential that we avoid this 
simplification. Internal consistency requires that we not 
impose such boundary conditions on the secondary struc- 
ture variables since we cannot impose them on the chain 
tangents in the bent configuration. Additionally, such 



periodic boundary conditions will substantially and ar- 
tificially reduce the statistical weight of the appearance 
of random coil (s = —1) segments along the chain. The 
cause of this artificial reduction is the following. With 
periodic boundary conditions domain walls must appear 
in pairs while in the physical problem they can appear in- 
dividually by destroying the native secondary structure 
from the ends of the chain. Periodic boundary condi- 
tions therefore suppress the probability of the creation 
of random-coil segments by one extra factor (3 (^ 1 in 
highly cooperative chains). 

The partition function is thus reduced to the remaining 
sum given by 



Z= (so\U(m)D N (m)U- 1 (m)\s N )r ls 



(11) 



where Dim) = U~ 1 (m)T(m)U(m) is the diagonalized 
transfer matrix and the remaining sums are over the sec- 
ondary structure of the first and last monomers of the 
chain as well as a single remaining angular momentum 
variable, m. We evaluate this final sums over sq,sn to 
find: 



i 



-(Af(m) + A£V)) 
Af(m)-Af(m) 



2^(i - Zm y + p 



(1 + /3z m - Z„ 



' h [l 



/?])]■ (12) 



Using the above partition function we may immedi- 
ately compute two measurable quantities: (i) the mean 
torque t(i/>) required to enforce the constraint on the 
chain tangents at either end, and (ii) the fraction of 
molecule in its non-native (random-coil) structure M(ip). 
The former can be directly measured by AFM in single 
molecule manipulation experiments while the later can be 
probed by circular dichroism spectroscopy. These quan- 
tities can be computed from the partition function shown 
in Eq. El by derivatives 



and 



r(V0 



M{i>) 



dhiZ{i>) 
dtp 



1 dlnZ{ij) 
"N dh ' 



(13) 



(14) 



Exploiting the symmetry u> m {n) = w_ m (re) we may col- 
lect the terms within the [■] in Eq. ^1 defining [•] = 
Z(m, h) to write the partition function as 



Z{ip) = Z(0, h)+2^T cos(ipm)Z(m, h) 



so that, using Eas. I13lf4*l we find 

—2 X)m=i 771 sin(f/'m)2(m, h) 



Z(0, ft,) + 2 YZ=i cos(^m)Z(m, h) 



(15) 



(16) 
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FIG. 2: a) The torque r(ip) required to hold the chain at a 
fixed angle ip for e w = 1.4 (solid line) and e w = 5.2 (dashed 
line) with h = 8. In both cases k> = 100, «< = 1 and N = 15. 
b) The corresponding number of segments having nonnativc 
(random coil) structure NM(ip). The buckling of the alpha- 
helix coincides with creation of one random coil segment that 
acts as a softer joint along the very rigid chain. 



and 



M(ip) 



dZ(0,h) 
dh 



+ 2E™ =1 cos(^m)^M 



N Z(0, h) + 2 Y% =1 cos(V>m)Z(m, h) ' 



(17) 



The series above are well approximated by partial sums; 
in practice, taking the first ten terms reduces the error 
to about one part in 10 6 . These numerically evaluated 
partial sums are shown in figure |2] At small values of 
the bending angle ip, there is a linear dependence of the 
constraining torque on ip. The alpha- helix bends like a 
flexible, elastic rod. At a certain critical angle ip*, how- 
ever, the constraining torque reaches a maximum and 
then drops precipitously for angles ip > ip* as shown in 
part (a) of figurc[21 This dramatic collapse of the chain's 
rigidity is akin to the buckling instability of a macro- 
scopic tube such as a drinking straw. The mode of the 
localized failure is, however, completely different. Ex- 
amining part(b) of figure |21 we see that at ip = ip*, M 
abruptly jumps to 0(1/N). The buckling of the alpha- 
helix is due to the creation of a single random coil seg- 
ment along the chain that provides a region of greatly re- 
duced bending stiffness. The size of the created random- 
coil section will remain on the order of Nk < / k > so for 
a large difference in bending moduli between the native 
and nonnative states of the chain, these "weak links" 
generically occupy a small fraction of the polymer. In 
the above example, there is only one weak link. 

The underlying cause of the buckling can be under- 
stood in terms of a comparison of the free energy in- 
crease per monomer associated with the creation of non- 
native states along the chain due to the reduction in 
chain bending energy due to the collapse of chain cur- 
vature into the more flexible random coil segment. At 



large enough imposed curvatures (at ip*) the reduction 
of chain bending energy more than offsets the free en- 
ergy cost associated with the creation of a chain segment 
having nonnative secondary structure and the buckling of 
the alpha-helix occurs. The mean conformation of chain 
changes from the arc of a circle (to minimize curvature 
energy) to two essentially straight segments connected by 
0(1) monomers in the random coil state where the ex- 
ternally imposed curvature localizes. It should be noted 
that critical angle ip* for the buckling of an alpha-helix 
depends not only on relative persistence lengths in the 
native and nonnative states of the chain, but also on the 
helix/coil parameters. If, for example, the chain coopera- 
tivity is decreased by changing e w = 5.2 — > 1.4 (dashed 
line to solid line in figure [2J) then the free energy cost 
of creating a weak joint in the chain is reduced and the 
buckling transition occurs at smaller imposed curvatures. 
One may observe similar effects through the change in h. 



The buckling transition may be considered as a type 
of nucleation process of random coil segments along the 
chain. Using this analogy, we can distinguish three dif- 
ferent types of buckling transitions based on the length 
of the chain and the chain cooperativity. For very long 
chains there is homogeneous nucleation in which the ran- 
dom coil segment first appears anywhere in the bulk of 
the chain, i.e. away from the chain ends. The creation of 
a random coil segment in the bulk of the chain requires 
the creation of two domain walls whereas the creation of 
random coil segments at the end of the chain requires 
only one. For long enough chains the added enthalpic 
cost of the additional domain wall is more than com- 
pensated by the translational entropy associated with 
the placement of that random coil segment. For shorter 
chains or chains having higher cooperativity, homoge- 
neous nucleation is replaced by heterogenous nucleation; 
the random coil segment appears at one end of the chain. 
Based on these considerations, the transition between 
homogeneous and heterogeneous nucleation of random 
coil segments should occur for chains where logiV ~ e w . 
Finally, even shorter chains or chains with still higher 
chain cooperativity, one can encounter a regime in which 
the entire chain spontaneously loses its native secondary 
structure at a critical angle. Such a transition should 
occur only if N < e w /h. The curves shown in figure |3 
correspond the intermediate case of heterogeneous nucle- 
ation of random coil segments. 



The appearance of the buckling transition is the first, 
and perhaps most dramatic consequence of the coupling 
the helix/coil, internal state variables of the chain to its 
conformational degrees of freedom. The generic conse- 
quence of such a coupling is the highly nonlinear bend- 
ing elasticity of the polymer as shown in figure |21 We 
now consider the effects secondary structure on the force 
extension relations of the polymer. 
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B. Stretching 

Before we develop the theory of stretching the 
HCWLC, we discuss the radius of gyration of such a poly- 
mer [43, |41j • To do so we note that the distance between 
the i th and j th monomers along the chain is given by 



3-1 



0*r 



(18) 



where j(s n ) is the length of a monomer measured along 
its mean chain tangent. Recalling that each monomer 
of the HCWLC represents enough amino acids (~ 3) to 
unambiguously assign a secondary structure to the seg- 
ment, the effective size of the monomer depends on that 
secondary structure. In the native, alpha helix state the 
monomer is shorter than in its generally more extended 
random coil configuration. To account for this aspect of 
the coarse-grained polymer model we define a monomer 
length that is a function of the secondary structure vari- 
able s n via 



7(s) 



7< 
7> 



if s = +1 
if s = -1 



(19) 



where, as the notation suggests, 7< < 7>, the length of 
a monomer increases when it loses its alpha-helical sec- 
ondary structure. We return to a discussion of reasonable 
numerical estimates of these values in the conclusions be- 
low. 



1. The radius of gyration 

We now compute the radius of gyration of the these 
polymers in solution by taking the thermal average of 
Rq N using Eq. ^| and the HCWLC Hamiltonian. Sep- 
arating the sum into terms diagonal and off-diagonal in 
the monomer indices we write 



N-l 

oat) = 

n=0 



N-l 



( 7 2 (s„)) (7(Sn)7(Sm)i 



n—0 m—1 



Both of these terms in the above expression can be evalu- 
ated in terms of the transfer matrices introduced in Eq.[S] 
We explore this point in some detail as it gives insight 
into the bulk of the calculations regarding stretching of 
the HCWLC. The reader who is uninterested in details 
of the calculation can pick up the discussion involving 
figure |3| in the last two paragraphs of this section. 

We note that the first (diagonal) term may be written 
in terms of the transfer matrices as 



4ir 2 Z 



e lm ^ N - eo) (s 0l m\D k {m)T 2 D N - k (m)\ SN ,m),(21) 



where both integrals range over the full unit circle. In 
Eq. |^ Z is the partition function defined in Eq. and 
we have also introduced the matrix L defined by 



L = 



7< 
7> 



(22) 



This matrix acting in space of secondary structure assigns 
the appropriate monomer length to the segment, i.e. 
(+l,m|r| + l,m') = (5 rn , )m /7< and (-l,m|r| -l,m') = 
8 m ,m'7> while both off-diagonal terms vanish. It repre- 
sents the action of the 'y(s) operator acting on a given 
state of the HCWLC, |s, m). The transition amplitude 
appearing in the above equation may be interpreted again 
as the amplitude for the fictitious quantum particle. In 
this case we compute the amplitude to propagate k imag- 
inary time slices at angular momentum m, be acted on by 
r 2 that measures the square of the length of k th segment 
and then propagate the remaining N — k imaginary time 
steps at the same angular momentum. Finally, the inte- 
grals over the initial and final angles of the chain may be 
performed explicitly; by not constraining these two end 
tangents we project out the in = state of the chain so 
that Eq. EH may be simplified to 



<7 2 (* fc )) 



^ z 

SQ,SN 



(s , m = 0\D(0) k T 2 D(0)"- k \s N , m = 0) 

(23) 



The remaining sums over the secondary structure of the 
initial and final chain monomers consists of only four 
terms and can be evaluated directly. 

Using the similar reasoning we may write the off- 
diagonal parts of Eq. [3UJ as 

(l(skhisj)h ■ ij) = 
2 E ^( S o,0| J D(0) fc r J D(l)-''- fe F J D(0) Ar ^| Sj v,0),(24) 

S ,S N 

where we have taken j > k. We observe that while the 
initial and final states are fixed at zero angular momen- 
tum (for the same reasons as in Ea. I23|l. the diagonal- 
ized transition matrix acting between the k th and j 
monomers is evaluated at an angular momentum of unity. 
To understand this we note that the thermal average of 
the scalar product in the above equation involves the av- 
erages of products of cosines of the form cob(0j) cos(0fc). 
These cosines generate m — ► m ± 1 transitions in the 
angular momentum basis so that the action of the cosine 
at the fc th monomer takes the initial m = state into ei- 
ther of m — ±1 states. The action of second cosine at the 
j monomer must return the angular momentum of the 
state to zero so that the integral over the initial and final 
angles does not cause this contribution to the transition 
amplitude to vanish. Thus of the four possible combina- 
tion of m — ► m±l acting at the two sites, only the two 
terms leading to no net change in angular momentum 
survive the final averaging. Because D{m) is even in m, 
the factor of two accounts for both of these terms. Once 
again the remaining sums over the secondary structure 
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FIG. 3: a) The radius of gyration of the HCWLC as a func- 
tion of the free energy cost per segment to transform to the 
random coil, nonnative state: h. In this curve «:> = 100, 
K< = 1, N = 10, and e w = 10. In the lower figure b) we plot 
the scattering function P(9) computed from the radius of gy- 
ration. The dashed line corresponds to that of a chain that 
is a mixture of helix and coils for h = 2.1 and e w = 10. The 
solid line and dashed-dotted lines are the scattering functions 
computed for the molecule in the all helix and all random-coil 
states respectively. 



of the initial and final chain segments can be performed 
directly. 

The final sums over j, k required to determine the ra- 
dius of gyration can be performed. If one were to assume 
translational invariance along the chain, the remaining 
sums over monomers in Eg. 1201 can be rewritten as a sin- 
gle sum 



N-l 



^2(N-k)C(k), 



(25) 



fc=0 



where we have denned the quantity C(fc) = 
(--f(si+k)j{si)ti+k ■ U) to be the tangent vector cor- 
relation function function weighted by the length of the 
chain segments. Due to our momentary assumption of 
translational invariance along the chain, this function 
is independent of the monomer index i. With this 
assumption the correlation function simplifies to 



1 



C(k) = -Tr 



r (oi) . D ^k . r (io) . £,( )JV-fc 
where we have introduced the matrices: 



r (or) 

p(10) 



U(0) 
U(l) 



-i. r . 
i. r . 



U(l) 
17(0) 



(26) 



(27) 
(28) 



written in terms of T defined in Eq. ^J] and the matrices 
U(m) that diagonalize the transfer matrix at angular mo- 
mentum to (see appendix IXjl . By performing the above 
trace, going to the limit of large N, and using fact that 
at any given angular momentum the helix phase is more 



probable than the random coil (i.e. Ai(to) > \2(rn) for 
all m), we may write the correlation function as the sum 
of two exponentials: 



C{k) = Ae- M « + Be 



(29) 



where the two lengths controlling the exponential de- 
cay of correlations are given in terms of the eigenvalues 
(Eq.DDJ} by 



log 



log 



Ai(0) 

Mi) 
MO) 
Mi) 



(30) 
(31) 



The coefficients of these two decays are written explicitly 
in the appendix IB1 

The existence of two exponential decay lengths for the 
correlation function reflects the fact that between any two 
tangents along the chain, the polymer may exist in one of 
two states having differing monomer lengths and thermal 
persistence lengths. To better understand this result it is 
instructive to study the limit in which it is highly unlikely 
to find the chain in the random coil state: e w ^> 1 and 
h^> 1. In that limit, A2(tt7) — ► so that If, — > oo; the 
decay of correlations is dominated by the length l a ~ 
1/k> in this limit. From extensions of this reasoning one 
can associate £ a with the decay length of correlations 
for a section of polymer that starts and ends with helical 
segments. Between these helical segments this correlation 
length is related to some function of both k<, re> due to 
fluctuations into the coil phase. Similarly, the length 
lb controls the decay of correlations for segments of the 
chain that begin and end with random coil segments. 
Once again, due to fluctuations into the helical phase 
this length is a function of both k < ,k > . 

Finally, given a form for the correlation function C(k) 
obtained from the assumption of translational invariance 
along the chain, we may directly evaluate the radius of 
gyration. We find: 



Ae 



«r) = JJ a _ 1)2 [JV(e*»-l) + . 



-N£ a 



1] + 



Be 



We suspect, however, that at least in the case of highly 
cooperative or short chains for which e w 3> log N the ef- 
fects of the chain ends will be significant and thus break 
the assumed translational invariance used above. In or- 
der to evaluate the radius of gyration of the chains that 
strongly break translational invariance, we numerically 
evaluate the requisite sums to determine the radius of 
gyration; these results are shown in figure [3] There we 
see the cross-over of a random coil (N > k<) having a 
Kuhn length of 7> = 3.0 to an essentially straight rod of 
length 7<./V in the helix phase. 

The radius of gyration can be experimentally probed 
via small angle elastic scattering. In figure 01(b) we plot 
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the scattering function for HCWLCs in dilute solution 
as a function of the scattering wavevector for scattering 
angle fi = ^ sin(£/2) In the same figure we show the 
predicted scattering from WLCs. 

A comparison of the predicted scattering HCWLCs 
and the better studied WLCs demonstrates that such 
scattering experiments alone are ineffective in differenti- 
ating between these two models. The scattering from 
a HCWLC can always be interpreted in terms of the 
scattering form a WLC having some effective persistence 
length. More generally, any measure of the radius of gy- 
ration will not distinguish the HCWLC from a simple 
WLC as long as the effective persistence length of the 
chain is adjusted to fit the data. In order to observe 
qualitatively novel behavior of the HCWLC, one must 
probe the force extension behavior of the chains. Here 
we will see highly nonlinear elasticity mirroring the non- 
linear bending elasticity of these polymers. 

2. Force extension relations: Small forces 

In the presence of a stretching force F the Hamiltonian 
of the HCWLC may be written as 

N 

H = H -FY / l(s l )cos(9 l ), (33) 

i=0 

where Hq represents the HCWLC Hamiltonian in the ab- 
sence of externally applied forces as shown in Eq. 0] A 
calculation of the chain partition function based on the 
above Hamiltonian would, of course, result in the com- 
plete description of the equilibrium force/extension re- 
lations for this model. Unfortunately, a closed form ex- 
pression for this partition function is not possible since 
Hq and the term proportional to the applied force are 
diagonalizable in the momentum and position represen- 
tations respectively. The basis states that diagonalizc 
the full Hamiltonian, Eg 1331 are the energy eigenstates 
of the quantum pendulum. We do not pursue this ap- 
proach here. The identical issue arose for the study of 
the stretched WLC; there approximate numerical diag- 
onalization j3|| and variational calculations [23, H|| 
have been successfully employed. 

We begin by considering small externally applied forces 
and consequently small chain extensions AL. Defining 
^7 = 7> — 7< to be the extension of a monomer under 
the helix-to-coil transition, we consider the small force 
to be those for which FA'y is small in comparison to 
the other four energy scales in the problem. Using this 
assumption, it becomes reasonable to expand the chain 
free energy in powers of the externally applied force. We 
thus generate a cumulant expansion 

oo 

logZ(V>)=£^ (34) 

1=0 

that is similar in spirit to those obtained from high tem- 
perature expansions of the Ising model |50| . In the above 



equation c/ is the I order cumulant. These cumulants 
are thermally averaged quantities in which the averaging 
is performed with respect to the zero applied force Hamil- 
tonian, Hq. Calculating derivatives of Eg . I3~H evaluated 
at F = allows one to calculate the mean extension of 
the chain in the direction of the applied force in powers 
of F. Finally, we note that since the remaining thermal 
averages are to be performed with respect to Hq, we may 
borrow the formalism used to compute the torque/angle 
curves and consider averages over restricted ensembles in 
which the first chain tangent is directed along F (taken 
to be in the i-direction) while the last chain tangent is 
fixed at an angle tp. We can thereby explore the coupling 
of applied torques to the extensional compliance of the 
chain using this formalism. 

Taking the first cumulant, which is the term linear in 
the force, we calculate the mean length of the polymer 
chain in the absence of any applied force 

N 

(L)(V0 = <I>(«fc) (35) 

k=0 

where the average (•),/, is taken over the force-free, re- 
stricted ensemble of chains having an initial tangent in 
the x direction and a final tangent making an angle 
■0 with respect to that initial tangent. The restriction 
placed on the first tangent breaks the rotational sym- 
metry of the system leading to a nonvanishing value of 
(L>(V0. 

Based on our discussion of the radius of gyration of 
the chain, we can compute these mean extension at zero 
force using the transfer matrix technique; Eo. l35l mav be 
written as 

N 

W = E E < s °' ° = 0|T fe 7 (sfc)cos^T w - fc | Sw ,^) w 

k=0 so,sjv 

(36) 

Working, once again, in the momentum representation 
we can recast the above expression into a simple sum over 
one angular momentum variable and the four possible 
combinations of secondary structure states of the initial 
and final chain segments: 

N-l 

(i)W= E E Vs N [(so,0\T k TT N - k \ SN ,l)cos^) + 

oo 

2 E (so,m\T k TT N - k \s N ,m + 1) cos([m + l]tp)]. 

(37) 

The structure of the above expression may be character- 
ized by using a simple graphical representation. In the 
left panel of figure 0] we represent the above terms for 
the first cumulant (the mean length of the chain) as the 
set of all one-step random walks in angular momentum 
space. As noted above in the computation of the radius of 
gyration, each factor of cosine increments or decrements 
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the angular momentum. In general the n th -order cumu- 
lant requires the determination of the thermal average 
of the product of n such cosines along the chain contour: 
(cos 9i 1 - cos 9i 2 ■ ■ ■ cos 0i„). There is a one-to-one mapping 
of such products to the set of all n-step random walks in 
momentum space. For example, the right hand panel of 
figure 0] shows all two-step random walks. That set of 
random walks in the momentum space sum to give the 
linear response of the mean length of the chain to the 
externally applied force F. 

By expanding about zero applied force we may obtain 
an expression for the series expansion of the mean length 
as a function of applied force of the form 



m + 2 



(L){F, VO = LoWO + LiWF + L 2 {^)F 2 



(38) 



where the prefactors of the odd- index terms, i.e. 
L2n+i{il>)F 2n+1 vanish upon averaging over all end an- 
gles, ip. Each term in the above expansion is the set of all 
n-step random walks in momentum space. Each involves 
a sum over the states of secondary structure at each end 
of the chain, which may be performed exactly and one 
infinite over a single angular momentum variable. The 
latter sum cannot in general be performed exactly, but, 
as discussed above, it may be numerically approximated 
to arbitrary precision. In practice because of the rapid 
convergence of this sum with high m, only a few terms 
are required to generate an excellent approximation. 

Constraining the initial and final chain tangents causes 
each walk in momentum space to be weighted by a phase 
factor exp(iAmip). Thermal averages over the ensemble 
of chains having unconstrained final tangents can be com- 
puted by averaging over ip. Due to the aforementioned 
^-dependent phase factor, this averaging eliminates all 
walks that results in a net change in the angular momen- 
tum of the chain. For instance the two walks comprising 
the mean length of the chain at zero force (left panel 
of figure both vanish when the final angle is uncon- 
strained. This is to be expected from basic symmetry 
considerations. By relaxing that constraint, one restores 
the rotational symmetry of the problem so that the mean 
extension of the polymer along the x axis necessarily van- 
ishes in the limit of zero applied force. 

Since the transfer matrix is even with respect to angu- 
lar momentum, there is an additional reflection symme- 
try; each walk from m to m! makes an identical contribu- 
tion to the final result as that walk reflected about m = 0, 
i.e. the walk from — m to —ml in which each increment 
of angular momentum is replaced by a decrement and vis 
versa. We employ this additional symmetry of the prob- 
lem to rewrite the set of two random walks comprising 
the mean length as shown in Eq. 1371 

We plot in figure the mean length of the chain as a 
function of h, the excess free energy per unit length asso- 
ciated with the existence of nonnative secondary struc- 
ture. If ft, < 1, we expect the chain to be driven into a 
random coil, nonnative structure in order to increase the 
chain conformational entropy associated with the disor- 
dering of the polymer backbone tangent vectors. For val- 



J k 

m + 1 



a j 



FIG. 4: Diagrammatic expansion used to calculate the ex- 
tension of the chain order-by-order in the externally applied 
force. The combination of two terms on the left give the 
O(F) terms that determine the mean length of the chain at 
zero applied force. The four terms on the right give the lin- 
ear response of the mean length due to an externally applied 
force. In both figures the angular momentum increases in the 
vertical direction and each horizontal leg of the walks repre- 
sent products of the transfer matrices at the labelled angular 
momentum. The labels 0, and N denote the beginning and 
end of each chain respectively. The labelled intermediate sites 
k and j denote the monomers where the cosines act to either 
increase of decrease the angular momentum of the walk by 
one unit. Mirror reflections of each diagram about the dotted 
line has the same contribution to the sum. 



ues of h > 1 we expect a highly cooperative (e w > log N 
in this example system) and therefore sharp transition 
to the alpha-helical, native state. This transition is evi- 
denced by the precipitous decrease of the chain's exten- 
sion occurring at h ~ 1. 

Using the formalism described above we may also com- 
pute the linear and nonlinear response of the chain to a 
force F by considering longer walks in momentum space. 
We report those results as a function of both F and ijj, 
the angle of final chain tangent. It is likely that in future 
single molecule force spectroscopy experiments it will be 
problematic to simultaneously control both the applied 
force and final chain tangent. While we suspect that force 
spectroscopy with unconstrained angles will be more ex- 
perimentally relevant, we believe that in order to discuss 
alpha-helical domain extensional elasticity within the na- 
tive state of a protein, such boundary condition prescrip- 
tions may prove necessary. We plot in figureEJthe nonlin- 
ear force/extension behavior (including terms up to F 2 ) 
of two representative HCWLCs for a variety of angular 
bends, ip. Higher order terms in the applied force can be 
computed similarly. For the case of unconstrained initial 
and final tangents, we have developed an automated pro- 
cedure to calculate terms of the perturbation expansion 
to arbitrary order |{5l|. 

In figure EJa) we observe the nonlinear extensional 
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FIG. 5: Mean length of the HCWLC in thermal equilibrium 
as a function of the parameter h - the free energy cost per 
segment to transform to the random coil, nonnative state. 
The bending modulus of the coil section k< = 10, k> = 100. 
and the length of the helix section 7< = 1 for this plot. The 
polymerization index of the chain is N = 10. The initial and 
final chain tangents are fixed so that 6o = &n = 0. As h 
increases we see the effect of the HC transition on the mean 
length of this stiff chain. 

compliance (through terms of order F 2 ) of the HCWLC 
for the case that chain is much shorter than its ther- 
mal persistence length. The different symbols correspond 
to differing imposed curvatures on the chain - see the 
caption. The finite extension at zero force is consistent 
with the assumption that the chains are simply bent into 
the arc of a circle consistent with uniform distribution 
of chain curvature as required to minimize chain bending 
energy. By noting the difference in slopes of the three ex- 
tension vs. force curves in figure Efa), one sees that the 
extensional compliance of the chain in the direction par- 
allel to the applied force depends on the angle of the final 
chain tangent, ip. In part (b) of this figure we observe 
the predicted force extension relations to the same order 
in the low-force perturbation theory for the HCWLCs 
where the persistence length is less than L even in the 
stiffer, alpha-helical phase. As evidenced by the equality 
of the slopes of all three figures, the effect of imposed 
curvature on the subsequent extensional compliance of 
this chain is minimal. For a more flexible chain twisting 
the final tangent has little effect on the compliance un- 
der subsequent extension. We expect the upper panel of 
figure to better reflect the mechanical of alpha- helical 
polypeptides. 

The perturbation theory presented above appears to 
be a useful approach to the study of small extensional 
deformations of alpha-helical protein domains. A generic 
feature of this model is the nonlinear growth of the ef- 
fective extensional compliance of the chain with increas- 
ing applied force. The underlying mechanism must be 
the enhanced statistical weight for finding a segment in 
its longer nonnative state. The applied force thereby 
accesses a reservoir of chain length built into each seg- 
ment; such a mechanism has been studied in detail by 
Tamashiro and Pincus 43] although that work ignored 
the role of the tangent vector fluctuations. These fluctu- 
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FIG. 6: The mean length vs. applied force for HCWLCs 
with constrained end tangents so that 8o — and 9m — ii. 
In both figures three different final angles are shown: -0 = 
(circles), tp = 7r/6 (triangles) and ifi = 7r/4 (squares). In a) 
n> = 100, k< = 1, e w = 10, h = 3.0, N = 10, 7> = 3, 7< = 1. 
In b) k> = 5,k< = l,e w = 10, ft = 1.5, N = 24. The up- 
per figure shows the expected behavior of an alpha-helical 
polypeptide chain. For comparison, a longer HCWLC of 
shorter persistence lengths is shown in b) 



ations should only enhance the nonlinear growth of the 
effective chain compliance. When k < <C k> the ran- 
dom coil sections of the chain recover a significant en- 
tropic contribution to their free energy. The chain as 
it extends into the random coil phase not only becomes 
physically longer, but also has a decreasing effective per- 
sistence length. The combination of these two factors 
enhance the effective compliance of the chain. 

The perturbatively evaluated HCWLC model, how- 
ever, fails to reproduce the central aspect of the high 
force limit. In order to capture this basic feature of the 
WLC that must also hold for the HCWLC, we must aug- 
ment our low-force perturbation theory by other meth- 
ods better adapted to the high force limit. Since one may 
compute the n th -cumulant, which generates the (n — 1) 
term of the series expansion of the force extension curve 
as shown in Eq. 1381 one might at first imagine that one 
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could perturbatively determine the force extension rela- 
tion of the HCWLC at arbitrarily large forces. This is 
not the case. In fact, the perturbative approach to the 
high force limit is incapable of capturing the essential as- 
pects of the behavior of the WLC or the more complex 
HCWLC. The high force limit, F — ► oo is an essen- 
tial singularity of the partition function and thus can- 
not be approximately by a Laurent series in force. This 
has been directly confirmed by Marko and Siggia j3|| 
who have shown by integrating out the transverse con- 
tour fluctuations of the WLC in the high force limit that 
the mean length approaches the maximal length, L as 
(L)/L ~ 1 — F~ x l 2 . Perturbation theory in the low-force 
limit (expansions in F) or in the high force limit (expan- 
sions in 1/F) will miss this result. 



3. High Force Limit: Mean Field Theory 

We now study the extensional compliance of the 
HCWLC in the high force limit. In order to explore the 
approach under high forces of the mean chain extension 
to its maximal extension L = 7> N, it is reasonable to as- 
sume that the tangent vector fluctuations become small 
so that a Gaussian approximation is justified. We may 
approximate the HCWLC Hamiltonian Eas. 141331 bv 

N-l , N 



N-l 



E - ^ - % }■ (39) 



i=0 



i=0 



The above Hamiltonian is now quadratic in the angles 
of the chain tangent vectors. We fix the initial and fi- 
nal chain tangents to lie along the direction of the ap- 
plied force, 9o = 9n = 0. By iterative Gaussian inte- 
grals over the remaining angles we determine an effective 
partition function that is now a sum over only the sec- 
ondary structure degrees of freedom. After integrating 
over 0i, ... , 6 N-i the partition function in the high force 
limit reduces to the form 

N-l 



TT E^ 1 ~ S k s k+l) + 



^ = E exp 

{«*} 

N 

■ J2(s k -1)-Fj2l(sk 

k=0 k=0 



k=0 
N 



J [{Si}} 



(40) 



where the remaining sum is over all 2 N configurations of 
the secondary structure variables. We have introduced 
the quantity J[{si},iV] produced by the Gaussian inte- 
grals. It is a function of that secondary structure config- 
uration defined by 



N 



J[{si},N] = H 



2tt 



2R 3 + k( Sj ) 



-7(so)F 



(41) 



Each term Ri in the above product is defined recursively 
by the equation 



FL — 



K(si_i)iij_i Fj(si) 



2R l _ 1 + K(*i_i) 



(42) 



for i = 2, . . . , N, where we fix the initial condition for the 
recursion by setting 



Ri = 



k(s ) + F 1 (s 1 ) 



(43) 



The i term in the product depends on the full set of 
secondary structure variables from site i — 1 back to 0. 
Similar recursion relations having a constant value of k 
and 7 are discussed by Lamura et al. [381 ]. 

Examining Eq. I4UI we see that by integrating out 
the tangent vector degrees of freedom we have taken 
the Ising-model partition function corresponding to the 
secondary structure variables, which had only nearest- 
neighbor couplings, and transformed it into the partition 
function for the secondary structure variables (sj) having 
interactions between the these variables at distant sites 
along the polymer chain. This result is to be expected: 
the combination of the coupling between the local chain 
tangents and secondary structure generated by k(s) com- 
bined with the long-range coupling of those chain tan- 
gents to each other over ~ k monomers leads to a new 
effective long-range interaction between secondary struc- 
ture variables mediated by the conformational degrees of 
freedom of the chain. It is clear from Eas. 1411 and 1431 
that the simple Ising description of the secondary struc- 
ture variables is recovered in limit of chains with a vanish- 
ing persistence length k — > where the tangent vector 
degrees of freedom do not mediate a long-range interac- 
tion between the Si variables. In that case the recursion 
relation can be trivially solved to yield: Ri = F-f(si)/2 
so that the remaining partition function of the secondary 
structure variables in Eq. 0U| reverts to that of an Ising 
model, but one for which each secondary structure con- 
figuration is weighted by its effect on the chain extension 
in the direction along the externally applied force. 

The short persistence length, decoupled limit is clearly 
not of primary interest in modelling an alpha-helical 
polypeptide. In fact, considering that we are primarily 
interested in molecules that are not significantly longer 
than their persistence length (in the alpha- helical phase), 
it appears physically reasonable to take a diametrically 
opposed approximation. For such chains where k > ~ N 
one suspects that the statistics of the secondary structure 
variables is better represented in a mean field approxi- 
mation enforced by the long-range interactions between 
these variables due to tangent vector correlations along 
the chain. To implement a mean field approach, we ig- 
nore boundary effects and study one secondary structure 
variable in the bulk of the chain, Sj. This single degree 
of freedom interacts with the mean field of all the other 
secondary structure variables along the chain. We define 



the mean value of these variables as m 



for all 
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j i. From this definition it is clear that — 1 < m < 1. 
In order to discuss the chain persistence length and effec- 
tive monomer length we must generalize Eas. l3ll9l respec- 
tively by introducing the mean values of these quantities 
by defining 



K = — (1 + m) + — (1 - m) 



7 



2 

7> 



(1 



2 

7< 



(1 



(44) 
(45) 



The linear dependence of these values on m may be jus- 
tified by noting that (sj) = m implies that each segment 
spends a fraction (m+l)/2 of the time in its native state. 
At least on time scales long compared to intcrconversion 
time between the Sj — ±1 states, one would observe the 
effective values R and 7 as defined above. Of course, 
nothing in the present analysis determines this intercon- 
version time, but we expect it to be on the time of con- 
formational changes of small molecules ~ 10~ 9 s. Both 
force spectroscopy measurements and protein conforma- 
tional changes occur on much longer time scales where 
the approximation Eas. 14414*51 is valid. 

We may write the mean field free energy of the chain 
under the externally applied force in the form 



^MF = -log] e-^^-V+i^- 1 ^^ 

Jmf [k, 7, m; N/2 - 1] • J[s h 1] • J MF [«, 1, m; N/2 - 1]}(46) 

We have defined quantity Jmf [k, 7, m; N] to be anal- 
ogous to J[{si},iV] (see Eq. fiOjl for a chain of N 
monomers with a fixed mean-field persistence length (k) 
and a fixed mean- field monomer length (7). We also take 
the mean field approximation: = m for all k in the HC 
part of the Hamiltonian. The function J[sj, 1] represents 
the one Gaussian integral associated with angular degree 
of freedom at the i th site. 

The physical meaning of Eq. is that the free en- 
ergy of the chain in the mean-field description may be 
written as the sum of three parts. The first part given 
by the negative logarithm of Jmf [k, 7, m; N/2 — 1] gives 
the free energy of the the half of the chain to the left 
of the selected site i. This free energy is evaluated using 
the mean field approximation for the secondary structure 
variables and the Gaussian (small angle) approximation 
for the tangent vectors. The last term in the sum is the 
analogous contribution to the free energy associated with 
the length of polymer to the right of the selected site i 
and evaluated using the same approximations. Finally, 
the middle term in the product appearing in Eq. 1461 is 
the contribution to the free energy of the i th site itself. 
The (Gaussian) integral J[si, 1] accounts for the tangent 
vector degree of freedom while the sum on one remaining 
St variable is written explicitly above. 

To justify choosing the i th site at the middle of the 
chain to be representative of any site, we must ignore 
boundary effects. Consequently the mean-field descrip- 
tion is most accurate in the limit of long chains, i.e. 



logN > e w . Finally, self-consistency requires that that 
the thermal average of the secondary structure at the i th 
also be equal to m. Thus we demand 



m = (sj) = 2 



MF 



dh 



(47) 



From Eq. 1471 we obtain a solution for m that, when used 
in conjunction with the mean-field free energy function, 
J"mf = — log Jmf [k, 7, m; N] gives a complete thermo- 
dynamic description of the chain under an externally ap- 
plied force. In particular we compute the mean length of 
the chain under these conditions from 



(L) = 



OF 



(48) 



We plot the mean extension of the polymer (L) in the 
direction of the applied force normalized by the maxi- 
mal extension of the chain L = Nj> in figure [7| The 
applied force has been nondimensionalized by the length 
7<. Qualitatively the figure may be discussed in terms 
of four regimes characterized by abrupt changes in the 
dependence of the mean length on applied force. For the 
smallest forces we see the initial extension of the pre- 
dominant native-state, alpha-helical chain. As long as 
the h > FA^ the free energy decrease associated with 
the breakdown of the native state of molecule and the 
consequent extension of each monomer by A7 = 7> — 7< 
is more than off-set by the free energy increase per unit 
length of creating segments in this nonnative state, h. 
Thus the secondary structure variables are frozen in the 
native state and the extension of the chain proceeds by 
the suppression of contour fluctuations transverse to the 
extension direction. The effective maximal extension of 
the chain is Lj < /j > (here L/3) and the saturation of 
chain extension reproduces the Marko, Siggia 36] result 
so that (L)/L - (1 - F- 1 /2)JV 7< / 7> . 

If the system were well-described by a WLC model, 
this high force plateau would be flat as the chain exten- 
sion asymptotically approached it maximal value. The 
observed slow growth of the chain extension or "pseudo- 
plateau" is due to the presence of local secondary struc- 
ture fluctuations. With increasing force, these fluctu- 
ations are biased toward the more-extended, nonnative 
state so that the mean length of a monomer grows slowly 
with applied force. This second regime characterized by 
the pseudoplateau terminates in a rapid extension regime 
where the secondary structure is pulled apart by the ap- 
plied force. The requisite force to open up these alpha- 
helices determines the transition to this rapid extension 
regime. That force is given by FAj ~ e w + h where 
the net extension of one segment enthalpically compen- 
sates for the creation of a random coil segment and a 
domain wall on the chain. Finally, in the fourth regime, 
the applied force has thoroughly destroyed the secondary 
structure. With increasing the force the now random-coil 
chain approaches its maximal extension L in a manner 
first discussed by Marko and Siggia. 
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FIG. 7: The mean length of the chain as a function of the 
applied force normalized by the maximum length L computed 
using the mean field approximation discussed in the text. The 
parameters for the solid line at e w = 10, and h = 1, and 
k> = 2 and k< = 1 and N = 10. The dashed line shows 
an analogous plot using parameter values of h = 1.5, e w = 8, 
k> = 100, and k< = 1. These latter values are representative 
of alpha-helical protein domains. 



IV. CONCLUSIONS 

We have proposed and explored an extension of the tra- 
ditional WLC polymer model in order to incorporate the 
presence of internal degrees of freedom along the polymer 
backbone and the coupling of those internal degrees of 
freedom to the conformational degrees of freedom of the 
chain. Such a model constitutes a minimal description 
of many biopolymer systems, but we focus on developing 
a description of an alpha-helical polypeptide chain. By 
studying the mechanical properties in thermal equilib- 
rium of such polymers one will develop insight into the 
mechanical properties of de novo designed alpha-helical 
chains and biomemtic synthetic polymers that acquire 
helical secondary structure in aqueous solution. 

Eventually one would like to apply such a model to en- 
tire proteins although such a program requires in general 
a mechanical description of beta-sheets and an investiga- 
tion of the mechanical interaction of protein subdomains. 
In principle, the combination of a model of the nonlinear 
elastic properties of such domains with an accurate three- 
dimensional protein structure in its native state should 
enable the investigation of protein dynamics and partic- 
ularly conformation change under biologically relevant 
conditions. It is known that at least some proteins are 
bistable having at least two structurally different confor- 
mations. In this paper we have shown that, due to the 
coupling of the conformational degrees of freedom to the 
internal, secondary structure variables, alpha helices are 
generically bistable mechanically. One may speculate as 
to whether this inherent bistability provides a mesoscopic 
mechanism to elucidate protein conformational change. 

In order to carry forward this program and to make 
quantitatively falsifiable predictions for the mechanical 
properties of single alpha helices it is necessary to deter- 
mine the for energy scales that enter the HCWLC Hamil- 



tonian. Unfortunately, these four energy scales are im- 
precisely known at best. The better studied energy scales 
involve the helix coil (HC) parameters e w and h. Both of 
these parameters are extremely difficult to estimate based 
on first principles since these energy scales involve com- 
plex solvation energies jE^ in addition to the formation of 
hydrogen bonds between adjacent turns of the alpha he- 
lix. These parameters can be estimated, however, by fit- 
ting HC models to both the results of molecular dynam- 
ics simulations and experiment |53l l54j . In terms of our 
parameters this work provides the following estimates: 
e w ~ 7 and h ~ 1.5. Thus we note that since typical 
alpha-helix domains in proteins have N ~ O(10) |55ll5^| 
these domains are highly cooperative, e w > \og(N). 

There is little data on the persistence length of alpha 
helices. We estimate the bending modulus of an alpha 
helix by assuming that its enhanced stiffness arises pri- 
marily from inter-loop hydrogen bonding. Taking the en- 
ergy scale of these hydrogen bonds to be ik^T - l^k^T 
|57| , an inter-loop distance of 0.36nm and helical radius 
~ O.lnm we find that (with k^T — 1) the persistence 
length of the alpha- helix should be i £ ~ 10 - 50nm. From 
this persistence length we determine the HCWLC bend- 
ing modulus via: ^p/7< = K > so that k> ~ 25 — 140. 
With the absence of hydrogen bonding in the nonnative, 
random coil state we assume that the persistence length 
is £p ~ lnm, typical of simple hydrocarbons [5^|. Thus 
it is reasonable to suppose that ip/j> — n< ~ 3. There 
is a significant dependence of thermal persistence length 
upon local secondary structure: the ratio k> / n < may 
be as large as 50. In order to explore the phenomenol- 
ogy of the model we have shown results for various pa- 
rameter values, but we have always included plots corre- 
sponding to these biologically relevant parameters men- 
tioned above. From these estimates of the biologically 
relevant energy scales in the model we predict the criti- 
cal torques for the buckling failure of the alpha helix to 
be ~ 40pN • nm. The force required to pull out the alpha 
helices leading to the dramatic lengthening of the chain 
is roughly 150 - 200pN. 

The central result of this paper is that an alpha-helical 
polypeptide is highly nonlinear in its response to applied 
forces and torques. The source of the nonlinearity is the 
coupling between local secondary structure and the con- 
formational state of the polypeptide backbone. Bending 
or pulling on the alpha helix mechanically can result in 
the abrupt breakdown of secondary structure and conse- 
quently a dramatic increase in bending and extensional 
compliance. The stresses required to access this highly 
nonlinear behavior occur on scales relevant to biological 
activity. 
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APPENDIX A: DIAGONALIZATION OF THE 
TRANSFER MATRIX 

The transfer matrix at a given angular momentum 
T(m) given in Eq. [S] is diagonalized by the similarity 
transformation: D(m) = U (pi) • T(m) ■ U(m) where 
the matrix U (m) is 



' \i(m)-2d m A 2 (in)-2d n 

U(m) = I 2c ™ 2c " 



1 



1 



(Al) 



where Ai i2 (m) are the eigenvalues of the transfer matrix 
given in Eq. 1101 and the functions c m and d m are simply 
the bottom row of the transfer matrix. In other words 



and 



r\ — P - h ~ K < T \k 

Urn c ± r n. rx, < 



(A2) 
(A3) 



APPENDIX B: CORRELATION FUNCTION 
COEFFICIENTS 



length squared, can be computed in terms of the matrices 
r (oi) j r (io) defined in Eqs. E3 EH1 by calculating the 
necessary trace, one finds that 



A = 
B = 



r (oi) r (xo) 
1 li 1 li 



r (oi) 

1 12 



,(10) 



(Bl) 
(B2) 



Determining these coefficients in terms of the fundamen- 
tal parameters of the model is now a matter of some 
algebra. To simplify this work and to better display the 
result, we find it helpful to write A and B in terms of 
7<.> and the transfer matrix eigenvalues Ai^ot) given 
by Eq.Hni We find 



.4 



4c c l7> - 7< (2rfi - Ai(l))(2dp - A 2 (0)) 
4c 2 l -(2d -Ai(0))(2do-A 2 (0)) 
4coc l7> - 7< (2d - A 1 (0))(2di - A 2 (l)) 



(B3) 



4c 2 



(2di-Ai(l))(2di-Aa(l)) 



and 



B = 



4ci7<(2rf - A 2 (0)) + c 7>(-2rfi + A 2 (l)) 
4 c 2-(2rf -A 1 (0))(2d -A 2 (0)) 
ci7<(-2d + Ai(0)) + c o7> (2d 1 - Ai(l)) 
4c 2 -(2d 1 -A 1 (l))(2d 1 -A 2 (l)) 



(B4) 



Here we discuss the coefficients A and B appearing 
in Eq. [2H These coefficients, which have dimensions of 



In the above equations we have used the functions defined 
in Eqs. IA2IA3I 
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